Phonon resonances in atomic currents through Bose-Fermi mixtures in optical lattices 
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We present an analysis of Bose-Fermi mixtures in optical lattices for the case where the lattice 
potential of the fermions is tilted and the bosons (in the superfluid phase) are described by Bo- 
goliubov phonons. It is shown that the Bogoliubov phonons enable hopping transitions between 
ferinionic Wannier-Stark states; these transitions are accompanied by energy dissipation into the 
superfluid and result in a net atomic current along the lattice. We derive a general expression for 
the drift velocity of the fermions and find that the dependence of the atomic current on the lattice 
tilt exhibits negative differential conductance and phonon resonances. Numerical simulations of the 
full dynamics of the system based on the time-evolving block decimation algorithm reveal that the 
phonon resonances should be observable under the conditions of a realistic measuring procedure. 
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I. INTRODUCTION 

One of the most intriguing prospects opened up by 
recent advances in atomic physics is the possibility of 
studying many-body quantum systems. In particular, ul- 
tracold atoms confined to optical lattice potentials have 
been shown to be perfectly suitable for implementing 
physical models of fundamental interest not only to the 
field of atomic physics but also to condensed matter 
physics [l|, |2|. Specific examples of a highly versatile 
many-body system include Bose-Fermi mixtures in opti- 
cal lattices, which have been used recently to analyze the 
effect of fermionic inipurities on the superfiuid to Mott- 
insulator transition [sHSJ. A further experimental setup 
closely related to condensed matter systems consists of 
ultracold atoms in tilted optical lattice potentials. Sev- 
eral fundamental quantum mechanical processes related 
to nonequilibrium transport of particles have been ob- 
served in this setup such as, e.g., Landau-Zener tunnel- 
ing 0, Bloch oscillations ^M, and processes analogous 
to photon-assisted tunneling [9|. 

Of particular importance, coUisionally induced trans- 
port of fermions confined to an optical lattice and cou- 
pled to an ultracold bosonic bath has been observed in 
an experimental setup of a similar type to the one con- 
sidered in this article [10]. Furthermore, a closely re- 
lated theoretical analysis of collision-induced atomic cur- 
rents along a tilted optical lattice was provided by Pono- 
marev et al. |ll| based on a random matrix approach. In 
essence it was pointed out in Ref. [ii| that the atomic 
current exhibits Ohmic and negative differential conduc- 
tance (NDC). 

In the present article we study a natural extension of 
the mentioned experiments, namely lattice Bose-Fermi 
mixtures with the fermions confined to a tilted optical 
potential. Motivated by earlier considerations [ij, [l3| we 
show that this system is highly appropriate for exploring 
the effects of electron-phonon interactions on nonequi- 
librium electric transport through a conductor. Specif- 
ically, we demonstrate that electron-phonon resonances. 



predicted to exist in solids nearly 40 years ago |l44l6l | , yet 
for which there seems to be no conclusive experimental 
evidence [i3i[l3> can be realized using ultracold atoms. 
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FIG. 1. (Color online) Fermions confined to a tilted potential 
V{x) form a Wannier-Stark ladder with a constant energy sep- 
aration A between adjacent lattice sites. The fermion-boson 
interaction enables phonon-assisted transitions from site i to 
j at the rate Wi^j. The released energy (i — j)/S. is dissipated 
in the form of a phonon emitted into the superfluid. 



In our model the fermions take the role of electrons, 
whereas the bosons provide a nearly perfect counterpart 
to acoustic phonons in solids. The tilt imposed on the 
lattice potential of the fermions corresponds to an ap- 
plied bias voltage to the system. The collisions between 
fermions and bosons result in fermion-phonon relaxation 
processes, which give rise to a net atomic current along 
the tilted lattice potential [Il|, [Ij, [T^ , as illustrated in 
Fig.[TJ We show that in our case, where both the fermions 
and bosons are trapped in an optical lattice, phonon res- 
onances in the atomic current occur. They arise as the 
momentum of the phonon emitted in the relaxation pro- 
cess approaches a so-called Van Hove singularity ^] at 
the upper edge of the phonon band. Moreover, we shall 
see that phonon resonances, at least on the level of ap- 
proximations made in this article, are also expected to 
occur in Bose-Bose mixtures I2li. 



There are crucial advantages of our cold atom imple- 
mentation over a solid-state system: First, neither impu- 
rities nor imperfections in the system suppress the res- 
onances in the current [17j and hence fermion-phonon 
scattering is the only relaxation process, which can be 
fully controlled via the bosonic system parameters p,|23- 
Second, large lattice tilts can be achieved with an energy 
mismatch between neighboring fermion sites that exceeds 
the bandwidth of both the fermions and the phonons. 
This makes it possible to study the influence of the pho- 
non density of states on the atomic current over the en- 
tire phonon band. In addition, it allows us to observe 
negative differential conductance [III [13, U^, which is 
realized in a solid-state system with difficulty by resort- 
ing to semiconductor superlattices [l^[l3,[2J|- Last, the 
parameters of our system can be chosen to ensure that 
Landau-Zener tunneling to higher fermion bands — often 
a significant effect in high field tr ansp ort [13 — is negli- 
gible despite the large lattice tilts [ij| justifying the use 
of a tight-binding framework. 

The key experimental techniques required for our 
scheme are twofold: first, the independent trapping of 
atoms of different species in species-specific optical lat- 
tice potentials |24l . |25| and, second, the tunability of the 
boson-boson and boson-fermion interactions by Feshbach 
resonances [^, [I^- Also, we exploit the possibility of im- 
plementing low-dimensional systems by strongly increas- 
ing the depth of the optical potentials along specific di- 
rections. In this way effectively one-dimensional systems 
can be realized by tightly confining atoms to an array of 
tubes [3. 

The structure of this article is as follows: In Sec. II 
we start from the Bose-Fermi Hubbard model and out- 
line the effective description of the Bose-Fermi mixture in 
terms of Bogoliubov phonons and Wannier-Stark states. 
In Sec. Ill we derive a general expression for the drift 
velocity of the fermions and show that this expression 
encompasses the phenomena of negative differential con- 
ductance and phonon resonances. Section IV contains 
the results of a near-exact numerical simulation of a re- 
alistic experimental procedure to determine the depen- 
dence of the atomic current on the lattice tilt. We end 
with the conclusions in Sec. V. 
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FIG. 2. (Color online) (a) The dispersion relation of acous- 
tic phonons in a solid with hcOq ~ lsin(ga/2)| (dashed 
line) and of Bogoliubov phonons in a superfluid confined 
to an optical lattice (full line) in units of the bandwidth 
w = 4Jt,-\/l + Ubnt/2Jt for the parameters UbTib/Jb = 0.5. 
(b) The corresponding density of states Dq — \dhu}q/d{qa)\~^ 
(full line) of the Bogoliubov phonons is approximately con- 
stant for small momenta and exhibits a Van Hove singularity 
near the edge of the first Brillouin zone (dashed line). 



where (J, j) denotes the sum over nearest neighbors. The 
operators a, {cij) create (annihilate) a boson and, sim- 



ilarly, the operators c 



create (annihilate) a spin- 



less fermion in a Wannier state localized at site j. The 
bosonic and fermionic hopping parameters are J(, and 
J/, respectively, and the on-site boson-boson and boson- 
fermion interactions are characterized by the energies 
Ub and Ubf, both positive and independently tunable. 
In contrast to the bosons, the lattice potential of the 
fermions is assumed to be tilted with an energy split- 
ting A > between adjacent sites. 

At low temperatures and for sufHciently small boson- 
boson interactions Ub/Jb ~ 1 most bosons are in the su- 
perfluid phase and accurately described by the phononic 
excitations of the superfluid. This description is obtained 
by transferring the bosonic part of the Hamiltonian Hbf 
into momentum space and adopting the Bogoliubov ap- 
proach [2^, which results in the phonon Hamiltonian 



Hph = ^ hujqb\bq 



(2) 



II. EFFECTIVE MODEL 

The specific system we consider consists of a homo- 
geneous, one-dimensional Bose-Fermi mixture of bosons 
and spin-polarized fermions, both trapped in separate op- 
tical lattice potentials. If the potentials are sufficiently 
deep so that only the lowest Bloch band is occupied then 
the Bose-Fermi mixture can be described by the Bose- 
Fermi Hubbard model [27] 



<i,j> 



<i,j> 



-C/b ^ «]«]«J«i + Ubf J2 "l"i4' 



(1) 



-J ' 



Here, q is the momentum running over the first Bril- 
louin zone, i.e., —n/a < q < n/a with a the lattice 
spacing, and the bosonic operators 6j (bg) create (an- 
nihilate) a Bogoliubov phonon. The excitation spectrum 
of the phonons is given by the Bogoliubov dispersion re- 
lation hujq — ^jEq[Eq + 2Ub'nb) , whcrc rif, is the bosonic 
occupation number and Eq — AJb sin {qa/2) is the dis- 
persion relation of noninteracting bosons in the optical 
lattice. The dispersion of the Bogoliubov phonons, al- 
beit not identical, is remarkably similar to the disper- 
sion of acoustic phonons in a solid-state system with 
hujq ^ I sin(ga/2)|, shown in Fig.[2la). The common fea- 
tures are a linear dispersion in the limit q — >■ and the 
band structure with a gap at the boundary of the first 
Brillouin zone. It should be pointed out that the accu- 
racy of the Bogoliubov description of phonons has been 



demonstrated experimentally \22\, thus we have perfect 
knowledge of and control over the phonons in our system. 
The fermions confined to the tilted potential are rep- 
resented by the eigenstates of the fermionic part of H^/ , 
i.e.. the Wannier-Stark states [22,,^]. The unitary trans- 
formation that relates the Wannier operators Cj to the 
Wannier-Stark operators dj is given by 



d,^Y.^^~^(^jf/^y 



(3) 



where J'niz) are Bessel functions of the first kind [31]. 
It follows immediately from the properties of the Bessel 
functions J^n{z) that the Wannier-Stark states are cen- 
tered at lattice site j and have a width of the order of 
Jf/A- By invoking the identities for the Bessel func- 
tions X]„ J7i-m{z)Jn-m'{z) = S,n,m' and {2n/z) Jn{z) = 

Sn+i{z) + Jn-\{z) one finds that the transformed 
fermionic part of iJ^y takes the diagonal form 



U^s = ^(jA)d]dj 



(4) 



The spectrum of il^s constitutes a so-called Wannier- 
Stark ladder [2^, [30| with a constant energy separation A 
between adjacent sites, shown in Fig. [1] We emphasize 
that the Wannier-Stark states are stationary and hence 
a net current of fermions along the tilted lattice only 
develops in presence of additional relaxation processes. 
In particular, it results from the inverse of the transfor- 
mation in Eq. ([3]) that a fermion initially localized at a 
single lattice site undergoes coherent Bloch oscillations 
with frequency A/S [3|. 

To rewrite the boson-fermion interaction in terms of 
Bogoliubov phonons and Wannier-Stark states we apply 
the same transformations as for the purely bosonic and 
fermionic parts of llhf- Invoking the results established 
in Refs. |28l . [2g | we find the fermion-phonon interaction 
Hamiltonian 



il^nt = UbfUb Y, d]d, + Y, [hA^''""' + h-c. 



with the corresponding matrix elements 



d^ d 
(5) 



fq,(. — Ubf 



f nbEq 

\NhUJg 



1/2 



^'(f^-?V 



i -iqai/2 



(6) 



where N is the number of lattice sites. We note that Hint 
coincides exactly with the description of the electron- 
phonon interactions in a solid-state system [29|| and thus 
is perfectly adequate for investigating related phonon ef- 
fects. 

The fermion-phonon interaction Hint describes the cre- 
ation of a Bogoliubov phonon out of the superfluid phase 
and the reverse process, both caused by the hopping pro- 
cess of the fermions. These incoherent processes involv- 
ing a single phonon enable fermions to make transitions 
between Wannier-Stark states separated by £ lattice sites. 
On the other hand, coherent processes resulting from 



Hint, in which a virtual phonon is emitted and reab- 
sorbed, lead to a phonon-mediated fermion-fermion inter- 
action, a renormalized fermion hopp ing and a mean-field 
energy shift of the fermions [12, 13] . 

The effective fermion-fermion interaction, being short 
range, can be safely neglected if we assume a filling 
factor of the fermions much lower than 1. However, 
we have to take into account the renormalization of 
the fermion hopping Jf when comparing our theoreti- 
cal model to experimental or numerical results. As ex- 
pected from theoretical considerations [l2l . Il3| and con- 
firmed by the numerical results in Sec. IV, the renormal- 
ization reduces the bare hopping. Last, the mean-field 
energy shift is exphcitly given by UbfUb + J^q l/g,oP/?i'^q, 
which is readily determined in the regime UbTib/ Jb ^ 1 
by taking the continuum limit p33 |. For a single delocal- 
ized fermion with J/ /A ^ 1 the shift is simply Ubf rib, 
whereas for a localized fermion with Jf /A ^ 1 we find 

Ubfnb{\ — aUbf/2hcs), with Cs = {a/h)y/2JbUbnb the 
speed of sound in the superfluid. 

It is essential for our model that neither interactions 
with the fermions nor the trapping potential confining the 
fermions destroy the phononic excitations of the super- 
fluid. These conditions are met by using species-specific 
optical lattice potentials [2J,|2a and by limiting the num- 
ber of fermions in the system l3J-l5j. Moreover, we restrict 
our analysis to fermions moving slower than superfluid 
critical velocity, which is Cs according to Landau's crite- 
rion [33] , to avoid excitations other than those caused by 
the hopping transitions. More precisely we consider the 
parameter regime C — 2aJf/hcs ^ 1, where 2aJf/h is 
the maximal group velocity of the fermions in the lowest 
Bloch band. 



III. ATOMIC CURRENTS 

Since the number of fermions is significantly smaller 
than number of bosons we effectively treat the superfluid 
as a phonon bath. Accordingly we describe the dynamics 
of the fermions by a master equation for the probabili- 
ties Pj that a fermion occupies a Wannier-Stark state at 
site j. The master equation reads 

dtP = Y [^^-(1 - P^WJ^^ - Pd^ - PJW^^J] , (7) 



where Wi^j are the rates for a phonon-assisted transition 
from site i to site j and the factors (1 — Pj) take the Pauli 
exclusion principle into account. The probabilities Pi ei- 
ther describe the occupation of a single fermion or the 
distribution of an ensemble of fermions; in both cases we 
choose '^j Pj = 1 . The average position x of the fermions 
at time t is accordingly x — '^Ao,j)Pj and the aver- 
age drift velocity Vd, which quantifies the atomic current 
along the lattice, is given by v^ = dtx = '^Aaj)dtPj- 

To obtain a useful expression for the drift velocity Vd 
we exploit the fact that the system is homogeneous so 
that the transition rates Wi^j only depend on the jump 
distance £ — i — j, where jumps with £ > are defined 



to be downward the tilted lattice, as depicted in Fig. [T] 
By using Fermi's golden rule based on the interaction 
Haniiltonian Hint, i-e., to second order in the coupling 
Ubf, we find 



Wf 



2tt 



^-J2\Ui\'iN, + l)S{Ju.,-iA), 



In 



w^ = ^Y.\^,A^N,&{^^,- 1^) 



(8) 



where N^ — (e' 



— (^^<iI^bT „ 2.)~i is the mean number of 
phonons with momentum q in the superfluid at tem- 
perature T and /cb is Boltzmann's constant. The rate 
Wf (M^/^) corresponds to the process, where a fermion 
jumps i sites down (up) the lattice and thereby emits (ab- 
sorbs) a single phonon of energy ^A. The expressions for 
the jump rates Wf and Wf" are valid as long as heating 
effects caused by the emitted phonons are negligible. 

To determine Vd in terms of transition rates we sub- 
stitute the expression in Eq. ([7]) for dtPj and find after 
some rearrangement that 
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FIG. 3. (Color online) The total drift velocity v^ (full line) 
as a function of the lattice tilt 5 = A/4J6. The contri- 
butions from different jump distances I (dashed lines) re- 
sult in a distinct NDC peak and a series of phonon reso- 
nances at the positions 5 Kj l/l. The system parameters are 
IJbTihl Jh ~ 0.5, Jf/Jb ~ 0.4 yielding for the phonon band- 
width 1.12 X 4:Jb and C, — 0.5. The drift velocity is given in 
units of vo — anbUbf/iJbh. 



Vd 



E 

j,i>0 



(ai) [<P,(1 - P,+e) - WtPj[Pj+i - Pj^t)] 



(9) 
with W^ = {2TT/n)Y,q\Iq,t?^'yTT-^q-^^) the phouou 
emission rate at zero temperature. The first term in the 
sum in Eq. (J9]) vanishes if the fermions are degenerate 
at zero temperature with all lattice sites j > jc occu- 
pied, i.e., Pj = 1 for j > jc and Pj = for j < jc with jc 
fixed by the number of fermions, whereas the second term 
represents finite temperature effects. For the system far 
from degeneracy and for small fermion occupation num- 
bers we may make the approximation PtPj ~ and use 
the condition ^ ■ P, = 1, which yields Vd = J2e>o(^^) ^e 
for the drift velocity. We note that at this level of ap- 
proximation the expression for the drift velocity and all 
subsequent analytical considerations apply, inutatis mu- 
tandis, also to lattice Bose-Bose mixtures |2l| . 

Taking the continuum limit of the sum over all phonon 
momenta (appearing in W^) and integrating over the first 
Brillouin zone we finally obtain 



a^N 



Vd 



E^ 



e>o 



dTlUJa 



\hA 



(10) 



q=qt 



with qi implicitly defined by Eq^ = ^^(fA)^ -|- {UhUhY — 
UhTib- The latter relation also yields the maximum jump 
distance £niaxj which is compatible with energy conser- 
vation in the fermion-phonon relaxation process. We see 
from Eq. (fTO| that the total drift velocity Vd is deter- 
mined by the sum over all admissible jump distances (. 
weighted by the phonon density of states \dTujjq/dq\^^ 
and the fermion-phonon coupling |/g/p. 

We first determine the mobility ji of the fermions in the 
Ohmic regime, which we define by the relation Vd — fJ-^ in 
the limit (5 — >■ 0, where S = A/4Jb is approximately the 
ratio of the lattice tilt and the phonon bandwidth for the 



small values of UbTib/Ji, we consider here. We find that 
Vd/S ex Y^Zi^^Je (^C) for (5^0, with C = 2a J f /he, 
as previously defined. This sum is a Kapteyn series of 
the second kind, for which a closed form is known [31j . 
yielding for the mobility 



M 



-U^f 






(11) 



The mobility diverges in the limit C — > 1, however, it is 
well defined in the regime C ^ 1- For the practically im- 
portant case, where the potential depth for the fermions 
and hence the hopping parameter Jf is varied, the mo- 
bility scales as /i ^ J? in the limit Jf ^ 0. 

The dependence of the drift velocity Vd given in 
Eq. (fTU|) on the full range of accessible lattice tilts 6 is 
shown in Fig. [3] We see that the atomic current makes a 
sharp transition from Ohmic conductance to negative dif- 
ferential conductance. The main reason is that the width 
of the Wannier-Stark states, proportional to J//A, and 
hence the overlap between states at different sites de- 
creases as the tilt A is increased. This reduced overlap 
results in a low drift velocity Vd — in analogy to NDC in 
semiconductor superlattices [1^ [l^ — and is determined 
by the matrix elements fqj. The drift velocity Vd de- 
pends on the tilt as (J//A)^^ for each jump distance £ in 
the limit of small fermionic bandwidth Jf/A <C f , which 
can be readily achieved in cold atom systems. As can 
be seen in Fig. [3] the NDC is particularly pronounced 
due to the superposed contributions from different jump 
distances i to the total current. These findings are in 
accordance with predictions for a free homogeneous su- 
perfluid [ill. Il3lll9| ■ however, the result for Vd in Eq. (fTU]) 
describes additional features stemming from the influence 
of the phonon density of states. 



Specifically, the drift velocity exhibits anomalies in 
the form of sharp peaks, which correspond to the an- 
ticipated electron-phonon resonances in a solid-state sys- 
tem [l4l.[la|. As shown in Fig.[2lb), the phonon density 
of states \d%bjq/dq\~'^ is approximately constant for small 
momenta, but exhibits a Van Hove singularity [20| at the 
edge of the first Brillouin zone. Therefore phonon res- 
onances arise as the momentum of the emitted phonons 
approaches the edge of the phonon band. It can be found 
from Eq. pTI)l that the resonance condition is given by 
(.5 = ^l + C/f,nft/2 Jb with (. = 1,2,..., f^ax and thus the 
current displays a peak at (5 ~ 1/ (. corresponding to each 
admissible jump distance L Since the fermion-phonon 
interaction provides the only relaxation process in the 
system these anomalies are directly reflected in the tilt 
dependence of the current. Further, the current vanishes 
as the lattice tilt exceeds the phonon bandwidth because 
there are no phonon states available in order to dissipate 
energy into the superfluid. 

Let us briefly discuss the effect of an additional (shal- 
low) harmonic trapping potential of the bosons on the 
phonon resonances. If the system is aligned along the 
X-axis the potential takes the form V{x) — mt,uj'^x'^/2 
with the trap frequency uj^. In the experimentally rel- 
evant case, where the harmonic oscillator length l^ — 
^JK/rnyjJx satisfies the condition Ix/N ^ nf,a, the lo- 
cal density approximation (LDA) is applicable |27l [34| . 
The LDA consists of replacing the bosonic density nf, 
by the position-dependent density ni,{x), determined by 
V{x)^ with the other parameters of the model unmodi- 
fied. Consequently, the average drift velocity in Eq. ([TU]) 
has to be replaced by v]f = j dxp{x)vd[nb{x)\, where 
p{x) denotes the normalized spatial distribution of the 
fermions. This averaging might cause some broadening 
of the phonon resonances. However, the position of the 
resonances 5 — i~^yjl + C/f,nf,/2Jf, and more generally Vd 
(except for the prefactor vq) depend on rib only through 
the expression UbUb/Jb ^ 1- Thus the trap- induced 
broadening can be made arbitrarily small by reducing 
the boson-boson interaction energy Ub, which sets the 
smallest energy scale close to resonance. This is partly 
explained by the fact that the resonances involve only 
phonons from the upper edge of the phonon band with 
wavelengths comparable to the lattice spacing a. 



IV. NUMERICAL SIMULATION 

The theoretical results derived in the previous section 
are, strictly speaking, valid for stationary atomic cur- 
rents in a homogeneous system. We now show based on 
numerical simulations that the predicted negative differ- 
ential conductance and the phonon resonances are ob- 
servable in a system of finite size under the conditions of 
a realistic measuring procedure. In order to measure the 
atomic current in an experimental setup we envisage a 
procedure consisting of the following three steps: 

(a) Initially, both the bosons and fermions are prepared 
in a horizontal optical lattice and the total system 
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FIG. 4. (Color online) The displacement of the fermion as 
a function of the lattice tilt 5 determined numerically (sym- 
bols + , o, x) and corresponding analytical fits (solid lines). 
The dominant NDC peak and the £ = 1 phonon resonance 
at the boundary of the phonon band can be clearly recog- 
nized. The ratio V between the height of the i — 1 resonance 
and the NDC peak is V = {0.09, 0.14, 0.18} for increasing J/. 
The £ — 2 resonance is visible for sufficiently large fermionic 
hopping Jf. The system parameters are J/ = {0.2,0.4,0.8}, 
Ub = 0.1, Ubf ~ 0.5 and the fitting parameters are e — 0.1, 
Jf — 0.7 J f, £max ~ 2, with energies in units of Jb- 



is in equilibrium [33 ■ The fermions are each local- 
ized in separate sites sparsely distributed through the 
fermionic lattice of sufficient depth so that PiPj « 
holds. In this configuration the fermions are auto- 
matically cooled by the surrounding superfluid p36| . 
which is only slightly distorted |13j . 

(b) Subsequently, the lattice of the fermions is tilted for 
a fixed evolution time t of the order of h/Jf to let 
them evolve. To obtain a detectable displacement of 
the fermions the fermionic hopping Jf may have to 
be increased by reducing the depth of the lattice. 

(c) Finally, the spatial distribution of the fermions is de- 
tected, e.g. , by in situ single-atom resolved imag- 
ing [33, [33 ■ From the difference between the ini- 
tial and final distributions it should be possible to 
extract a reliable estimate for the atomic current. 
Alternatively, the momentum distribution and hence 
the drift velocity of the fermions may be determined 
directly by a time-of- flight measurement. 

To demonstrate the feasibility of this procedure we 
have simulated the fully coherent dynamics at zero tem- 
perature of both the bosons and fermions based on the 
complete Bose-Fermi Hubbard model given in Eq. p]). 
To this end we have used the time-evolving block dec- 
imation (TEED) algorithm [3^, which is essentially an 
extension to the well-established density matrix renor- 
mahzation group (DMRG) method [13, Ell. The TEBD 
algorithm permits the near-exact dynamical simulation of 
quantum many-body systems far from equilibrium, which 
is crucial for our purposes. The experimental procedure 
was simulated for a system consisting of f Of lattice sites 



with a bosonic filling factor rib — 1 and a single fermion 
initially located at the center of the lattice. We used box 
boundary conditions and the evolution time was limited 
to T = 24?i/Jf, to minimize finite-size effects, primarily 
reflections of phonons from the boundaries. 

Figure |4] shows the displacement of the fermion as a 
function of the lattice tilt 5 for a set of realistic experi- 
mental parameters, in particular, for three different val- 
ues of the fermion hopping Jf. The main features pre- 
dicted by our theoretical analysis can be clearly recog- 
nized, namely the dominant NCD-peak, the phonon res- 
onance at i5 ~ 1 corresponding to the jump distance (. = 1 
and the suppression of the current for (5 3> 1. In addition, 
for sufficiently large values of Jf the phonon resonance 
at (5 « 1/2 corresponding to the jump distance i. — 2 \& 
visible. For lower values of Jf the phonon resonance for 
1 = 2 may still lead to a characteristic drop in the current 
once the lattice tilt exceeds half the phonon bandwidth. 
The resonances for higher values of the jump distance i 
seem to be masked mainly by broadening effects for the 
parameter regimes tested. 

The ratio between the height of the £ = 1 pho- 
non resonance and the NDC peak, i.e., their relative 
visibility V , depends nontrivially on the fermion hop- 
ping Jf. The hopping Jf enters the expression for Vd in 
Eq. p^ through the matrix elements /q_^ in the form 
of J^i{siJf/Ji,), where the parameter s depends on the 
lattice tilt. Close to resonance we have s ss 1, whereas 
for small tilts £A ^ Ubn\j we obtain s « yj2J\JU\inh and 
hence s ^ 1 in the superfluid regime. Accordingly, the 
height of the ^ = 1 phonon resonance (s ~ 1) varies only 
slowly with J f as opposed to the height of the NDC peak 
(s ^ 1), which is characterized by the oscillatory nature 
of the Bessel functions. Thus a careful choice of J f al- 
lows us to optimize the visibility V ^ i.e., to minimize the 
height of the NDC peak, by tuning siJf/Jb close to a 
zero of the relevant Bessel functions. For the increasing 
values of Jf considered in our simulation the first zero of 
the Bessel function J^i.'z) is approached, which explains 
the improved visibility V for higher J/, shown in Fig. |4l 

The general broadening of the phonon resonances may 
be caused by finite-size effects and multiphonon pro- 
cesses. The finite-size effects include reflections of pho- 
nons from the boundary of the system that in turn af- 
fect the motion of the fermion. Multi-phonon processes, 
which have been neglected in our theoretical analysis, 
also cause some broadening of the single-phonon reso- 
nances. In particular the continuous drop in the current 
to zero can be explained by the emission of several low- 
energy phonons during a jump process, which is allowed 
even if the lattice tilt exceeds the phonon bandwidth. 

In order to compare the numerical and the theoreti- 
cal results, at least qualitatively, we adapt our expres- 
sion for the drift velocity Vd. to a system of finite size. 
In addition, we introduce the observed broadening ef- 
fects characterized by the energy e into the theory. Ex- 
plicitly, we calculate the average drift velocity by us- 
ing the expression w^ = 'Ylit>^ ('^^) ^/ with the rates 
Wf = (2^/V^£;i)E, \UA^ exp[-(?ic^, - i^fle\ The 
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FIG. 5. (Color online) The displacement of the fermion as 
a function of the lattice tilt 5 for different evolution times r 
determined numerically (connected data points). The NDC 
peak and the 1=1 phonon resonance are clearly visible also 
in presence of strong boson-ferniion and boson-boson inter- 
actions. The additional peaks between the NDC peak and 
the phonon resonance are due to finite-time effects. The dis- 
placement increases approximately linearly with the evolution 
time r. The system parameters are J/ = 0.5, (/(, = Uhf = 1, 
with energies in units of J;,. 



rates reduce to the previous expression W^ in absence 
of broadening, i.e., Yan^^i^Wf = W^*. The average dis- 
placement of the fernrion (in units of lattice sites) af- 
ter the evolution tinre t is then approxinrately given by 

We fit this theoretical model to the numerical results 
with the broadening e, the maximal jump distance £niax, 
and the effective fermionic hopping Jf as free param- 
eters. The first two parameters allow us to extract a 
quantitative value for the broadening and to determine 
the donrinant hopping process in the experiment. The 
fernrionic hopping needs adjustment because the bare 
hopping Jf is renormalized by coherent phonon pro- 
cesses, which lead to a reduced hopping Jf < Jf as dis- 
cussed in Refs. (l2l . [l3| . This is in direct analogy to the 
increased effective mass of polarons due to the drag of 
the phonon cloud |42| . 

As can be seen in Fig. |4]the fit describes the numerical 
results very accurately with moderate broadening e and a 
minor reduction of the fernrionic hopping Jf. Further, we 
find that Imax = 2 and thus the donrinant transport pro- 
cesses are nearest- and next-nearest-neighbor hopping; 
this is consistent with the absence of higher-order phonon 
resonances noted earlier. The high level of agreement be- 
tween our theoretical nrodel and the nunrerical results is 
partly explained by two observations: First, the numeri- 
cal results show that the fermion reaches a constant drift 
velocity on a time scale much shorter than the evolution 
time T, thus transient effects due to the sudden tilting of 
the lattice are negligible. Second, the drift velocity re- 
mains approximately constant over the entire evolution 
time, as illustrated in Fig. [S] resulting in a stationary 
atomic current. 



Since our numerics simulates the full dynamics of 
the system we can monitor the average position of the 
fermion at different evolution times r and extend the 
range of parameters considered so far. In particular, we 
are interested in the atomic current in presence of strong 
boson-fermion and boson-boson interactions, where our 
theory, assuming a superfluid phase, is no longer applica- 
ble. Figure [5] shows the displacement of the fermion as a 
function of the lattice tilt 6 for the interaction strengths 
Ub = Ubf = Jb and for different evolution times r. We see 
that the NDC peak and the £ = 1 phonon resonance are 
not noticeably affected by the presence of strong interac- 
tions, thus they both seem to be robust features of the 
atomic current. We note that the additional peaks in the 
displacement-tilt dependence (located between the NDC 
peak and the £ = 1 phonon resonance) are due to finite- 
time effects, which is revealed by a more detailed analysis 
of the time dependence of the numerical results. Further- 
more, we observe an approximately linear increase of the 
displacement of the fermion with the evolution time — a 
finding also confirmed for weaker interactions. 



Stark states, respectively, with a generic fermion-phonon 
interaction of an identical type to the one encoutered in 
solids. 

We have studied the dependence of the atomic cur- 
rent on the lattice tilt from first principles and found 
that our model accommodates negative differential con- 
ductance and phonon resonances. To demonstrate that 
these features are observable by using ultracold atoms in 
the context of a finite-size system and a realistic mea- 
suring procedure we have calculated the atomic current 
numerically by using the TEBD algorithm including the 
full dynamics of both the bosons and the fermions. Our 
numerical results show that the phonon resonance at the 
boundary of the phonon band is a robust phenomenon 
that occurs over a wide range of system parameters de- 
spite broadenin g, w hich might be increased by finite tem- 
perature effects [13, Il3 not directly taken into account in 
our model. Finally, we note that in an obvious extension 
of this work we will investigate the effect of the transition 
of the bosons from the superfluid to the Mott insulator 
regime on the atomic current through the system. 



CONCLUSIONS 
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